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(54) Doppler radar 

(57) There id provided radar apparatus, and a meth- 
od for determining the approach velocity of a target 
whereby the target is first classified as a "target of inter- 
est- by a thresholding process applied to successive 
signals to identify potenM^^ 
is then applied only to aselected portion of data relating 
to the target d interest. 

This invention enables existing moving target mA- 
cator MTI radars (1)to be upgraded to give full Doppler 
capability wKhoutthe excessive processing capacity re- 
quired by conventional Doppte' radars - 
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Description 
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Tt» pm*t invpnta -.fetes to otmm '« wd » ™"»" «• *» DSi'lZL 

the time token by the beam (of finite width) to pass across the target position. This limits the ™ m "™**™ 

from all strikes constituting one look at a target must be m some way cohe, entty Integrated when measure ^ 
shift of the target in question. aM1 ,^ 1 4 i ,,^titf irtsteadthev are detected by obsenring 

imppBIre a**, open. IhP possibility =1 •»*»» m-M"™"" ^^7*^ "s2 ?£ HS™. A THrcsurod 

^SpT^=s^ 

p U tee basis. This is known as PRI (or PRF) stagger (PRF = P"* 8 ™^^^^ v 
emptoys a short bu*t of pulses at one PRF followed by a second ^^f^J^^^ of two or 
arTtransmitted altematefy in a regular manner. Strike returns from a ^ J« burst are 

"ore different bursts in order * ensurej 

anaVsedtodetenmme the apparent pnawchange M ^ J^V™^™^^ anaMttocirtwininaQB^Mhlc*! wm. 

correct Doppler shift and( . |imited bv dLltter echoes and noise, by the (relatively email) 

in practice, the measurement accuracy of e, and ^ is hrmteo oy ctiraer «=™ ^ 

numbed o. strikes per burst, and by stabilKy -P^^' n *^ '^S^Z^**™ match — 
calculated from * and ^ are therefore eubiectto smaU e J^JJ*^^,^ ^^d. However, 

butrnore.hanoneve.uerr^ 

bycarefuleeteouonof the - drtda»P»»»^2 , "'V fc range tor targets ol interest. 
aTe well spaced in .regency, so that only the correct value he ^ J^SSl-^L* enough details of im- 
«• Similar principles apply to radars ^f^^^'^^Z^t^^)^ Dopier processing 

pi.montaticn are. of course, different. In both ^^^^^^1^ into different Doppler bins 

^Sc^ 
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> w*rw Pffactit/e and nol oniy positively identities a moving target from • 

c ,uu:^ 

each of which requires a complete ^'P^^^^dteBnqgw movingto^te-charaoterteedbynon- 
One very important advantage of Doppler radars ™ * 3Uch as toU buildings, terrain features, 

zero Dopptershlfts, Jromsignats^ftected from ^^ a ^l^^^^^Q Wg6tB ^ ^i are detected 
XS.etc.fa^gwKhinlhe^ are aUocate d. the clutter 

in differ** Doppler bins to those to whjch ^^f^^S^SZ^ ***** » * U8ed tQr rapidly 
echoesdo not compete with ^2Xit^^d disregard theclutter 

The essential principle of MTI processing ©to ^*|^V** succ^s'we strifes on stationary clutter or slow-moving 
PHI is likely to be. at most, a tew millisecs. e ^^^^^^ nV a smell residue. By contrast, sue- 
targets should be almost identical. ^^^^^^^^^ received. Subtracting these 
cesstve strikes on a rapidly moving target wM show a appfoachBS K (i BO degrees) the 'sub- 

strikes returns should leave a much larger residue «^ r ^ P ^%^^m a n either signal on its own, 

^owover.it.e.ppr^ve^o.amov^ 

pulses is exactly 2K(or±2N* etc) Hs echo be f™^^^^ radar. 

will be lost. This is termed Doppter blmdness and » cteery _analogou ^ ^ » ftjthcr pRp bufSt 

The standard cure for Doppler blindness « »^ P"^^ J^SSd be seen ctearty when a dHlerant 
mode transmissions, as in a Doppler radar: ^T^^^SSnV^ <*°* 3 d <*" ^ 
PRl is used. The method is not perfect but as with ^SbVe for ^geis of interest, 

the towest residua, biind speed can be raise to f^^^S^Shat of a Doppier radar having 

in summary, an MTI radar can achieve clutter ^J^^^S^. stabilrty and PRF stagger attributes 

similar detection performance. The MRI ladar ^^^^^XimXed and. if desired, can easily be impte- , 
slmifcrtomoseaatfcpp^ 

mented without using advanced digital techniques. M 1 1 * eon "° 1 ?^* te 

reliability at modest cost; very many ™<* J^*?'* requires less processing power than present 

taf9 tcorang to the flrst aspect * tha Invention there is provided radar apparatus comprising: 
a traiismmer arranged to transmn a serfes of^lses; 

a receiver arranged to receive and store successive ^ si ^ ™^ process applied to data originating 
I Sector arrayed to identify ^fJ^^^^^^SSxo process stored data 

potential target corresponding to that echo. 

By employing the present invention it b ^^^^^^^^r^dS^m to the target 
or target interest and then, once such a ^^^^^^Lto^6^ by the Doppter 
in order to obtain the approach velocHy. ^^S" ™^ target which wHI be distinguish^ 
extraction processor. An "echo of mte.est ^^^J*"*^^ ^ e , J,ake Wo account other parameters 



oeoverwrmen.u»Biei"'=«"'7 ■»—■■■" — _ . , 

chronologically which may be achieved by t.me taggmg the data. 
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In this manner the Doppler extraction processor can have access to the data required to assign a« approach 
velocity to a target of in.eresi. without having to have excessive memory storage. The data can conveniently be stored 
in RAM in the Jorm of complex video comprising (p) and (q) signals, one « Phase andone in *>*^ 

Preferably the detector comprises means tor subtracting successive ,ece.ved echoes andthret folding the drffe 
s ence signal obtained to determine an echo of interest. Thresholding processes can be implemented us.ng relatively 
little processing power compared to periorming Doppler processing every point on an «™aa 

Advantaoeouslvthe Doppler extraction processor Is arranged to: retrieve from the stored data the relativephases 
of s^S^^omTurget of intere* compare a ran* of velocity vah.es with the 
Senile velocity value that belt fits the phase data. (The range* velocity values .sprefer^^ 
.c ranoeta -a eartJlar target type.) The be£ fit can be found by calculating a fit error value for each velocity value n 
J^VSSSSS^X vaiues constituting an error function; prominent nulte in '^function arable 
values of the target approach velocity. The approach velocity correspond^ to the deepest null ™ th £*»j"f£^ 
toa teraet of interest provided the null depth is below a predetermined threshold. Tocompare null deptf* more accu- 
rateiSeDCler^ 

steptore^g^^ 

S the t unX In the vicinJy of these nulls, using a relatively fine step hi order to determine which c the coarse nulls 
is in tact me deepest and whether it crosses the final acceptance threshold ^tahted -least sum 

so in an altemaUve to the above, the Ooppter extraction processor can be arranged to apply a wealed toast sum 
of Z£* TZ nuN to determine the nXmum of the null. Performing a weighted ' «" u ^ 

£^ rioortthrrvrt enables the absolute minimum nuO value to be determined for each candidate null regwn. One 
2^^Xi^X^' algorithm b not based on any assumption that the ref erence phase * « of t» 
22X11 has beeTmeasured with negligible error. This allows deeper and more rehab* null values to bo found 

2s in cases where the signal to noise ratio (SNB) Is relatively poor. . . k „„ t nll i ftaa 

*t is preferable that the pulse repetition interval is varied, efther on e pulse by pulse bas.s pr on a burst of pulses 

basis? since this reduces ambiguity in the Doppler values. taroat 
oaa .s, rv->n D | S r extraction processor may advantaggously be used to verrfy the target as txans a "target 

of *I2Zffi2212£2222^^ ** vd,oc,,y value teeK 10 * ou,put bu l 

*J222 to cSs ^oweve? Doppler measurements obtained according to the Invention can a^.beexptorted 

22!!2!*&21^T* 

Failure to extract a clear Doppler measurement for a reasonably <^ detect cn not suppressed by MTI thus ndicate* 
that the taroet is orobablv due to distributed moving clutter, ie. either false or unwantea. .^i^^, 

by an MTI thresholding process. ^tk*rote rt ™id^ar^ 
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too*™ » b. «*« r™. «. M» p™— "> — « OT 
of a target, the method comprising: 

Th. m. «., » wayof —«»•* -.——«.•>. «™^«g»re..<-«* 

Figure 2 auwa»s«n«i™pl»yp»«»™p ,tt ' 0nln,e ' v * l: . . 

F ^.3«u^«»l^^9»*^ a,M " ma ™ l,m ' d ** ,i ' S °'°* 0 " : " < ' 
F*». 4 M». WW pM»«n»19»«y *» » ""<""«' « M " < * **" 

H» Invention is described b»toW » ^^"^i^ ^^g^^Jr^^^U^S^toO^Mm^t^^'^"^^^^^ 

management algorithms. a treauencv synthesiser and synchronising unit 3. Because the 

InRgure t the existing MTl ^ 
MTl function requires the radar transmitter 4 and nece'ver 5 » measijr H emB nt system makes use of PRF and 

pulses an MTl radar will Irifclude a synchronisation una 3 . The Doppler measurem 

other synchronisation signals generated by theunita d| d unh 6 . Efficient clutter can- 

The receiver S supplies echo senate m J^^^eSeTdouble or triple cancellation per- 

cellation. avoiding both Doppler blindness ^ ^ase bLr^ rc^*W . ft and me ^ (q) h quad- 
formed on complex video s^Ccm^x ^SCSKI^^^SLl la some mannet However, 
rature. In some MTl radars the two ^^^^^SSSSL signal are presented in separate channels. 

pulse compression is used. MH «™*d b« either analogue or digital means. Digital implementation is nor- 

MTl subtraction processing may be performed ^J^^^L^ ^ |geo 6ignate . The invention will therefore 
many preferred where multiple subtraction » to be " SSiSSa Signal samples rate of. say. 10 MHz 

-:^rrrrdis 

synchronisation unit 3. . sitinals rt ua to stationary clutter and slownrnoving targets 

The MTl canceller and detectoruml 9 J^^^S^Si^ <* '» M to * '* mU8 ' 

before thresholding the residual video to obta,n d ^^°? m ^ a meanB there is a short time lapse between 
and subtract video echoes received over -several ^l^f^^ nat ^^^ci^W^^^ 
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by parameters (e.g. azimuth, range bin. amplitude) expressed in unite directly related to reference datum* (toning 
pulses, etc.) used wrthin the detection process but not necessarily meaningful to the end user. 
P The raw detection data passes to a plot extractor unit 7 . This unit stores and associates <T« 
algorithms which deckie whether adjacent detections shouW be reported ^^^^^^^!St 
detection report Unreliable, implausible or uninteresting detects may be rejected and the remaining target s of 
bVreportec^ ttnjt plots in a data format and related to a coordinate system which ,s meamngful to the 
end^er The £Z Z» £r £-1*. be a dispfcy console, a target tacRing processor, or a muft-sensor date 

iUSi Z STperform Its tasks efficiently the plot extractor unit must temporary store raw detection date and win 
*J2£S2£ some latency additional to that introduced by unit 6. A feature of the ivenfon . that ,t exploits the 
latencies inherent in the existing detection and plot extraction systems. . iiri . . _. rlltp . „ 

The invention uses Doppler measurement algorithms which process a sequence of echoes returned 
the tar 9e rrce^plel7hlft is being measured Since these echoes arc presented sequentially in the vdeo sgnal 
*L SroSL3oxinnaiely equal to the radar dwell time, ills necessary to store the video signal temporanly. for 

teSSS^^SK^wg-. Video signa* are already stored tennporarily wthin -unit * 
tenteSo^ 

unit it is preferable to provide an additional cyclic video store 8 for Doppier processing. MMabutea of 

si^f the dtettal vkteo sampling rate is typically around 10 MHz this store requires at most a few Megabytes of 
<JZ^^E£ZXS?L* c-Lucus-y written into store, signals whfch have beenin store formore 

hoori 11<;pd w Doooler processing. Tne Doppter processes management unit 9 receives uming 

quantitatively where this is helpful and to convey some estimate ol typical values, the description g 
to a fictitious radar having suitably rounded but fictitious operating parameter values. 

For a nominal radar frequency of 1000 MHz (= 1 GHz). 
The nominal wavelength is (X) approximately 30 cm (= 0.3 m). 
30 The standard formula for Doppier ahHt F(d) on a moving target echo is. 

F(d) = 2xV/(X) 

ss wheraVismeradiaiveloc^ 
th9 radar). 

It V = 150 m/sec (typical of a subsonic aircraft); 

F(d) = 2 x 1 SQ/D.3 = (+)1000 Hz. 

Note that the Doppter shift is very small indeed by comparison with the radar fr •gj^g^jj* ^ ***** 
Z or D^plefn^suremen, f-«^2^^SSKS5 ESSES Z£l ^ would 

is: 

1/0.03 = 33 3 Hz; this corresponds to F(d) = 5 m/sec at A. = 0.3m. 

Adap^theradarforD^ 
■ocity bSmfeee. This degree of resolution ^^^SS^ri^lin^ 
for example, Btow-movlng ships from nea*y eeaclutt er. J^Wng echoes from slowing 

processing load whilst maintaining a required velocity resolution. 
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Consider a sequence ol strike echoes received during 30 msec* dwell time on one look al an isolated target. 
(Cases where targeis ere not totaled will be considered later.) During th* dwell time the target will move d.5 m closer 

to the radar by virtue of its approach velocity V= 150 rrVsec. . ^ . ^ . Martiti ^ ani 

If the radar transmits 5 MHz pulse compression bandwidth and compresses the received echoes to an equivalent 
5 pulse width of: 

compressed pulse width (cpw) = i/(bandwidth) =0-2 rnicrosecs. 
to The nominal range resolution wiD be: 

range resolution = (cpw) x cV2 = 0.2 x 30072 = 3Qm [c = 300 x 10*6 m/sec]. 



iS 



20 



26 



Thus the effective width ot each range bin used in the radar receiver cannot be less than 30 m; this is many times 
areater than the distance (4.5 m) which the target will move during the dwell time. 
9 icanme^ 

useful echo signals from the remaining strikes ot the sequence. Moreover, this assumptlcn w,.L m gen e ral. be valid tor 
taroat velocities ud to V =±1000 m/sec, Le. for all targets likely lobe ot Interest. 

^^^^r^ processing only those compressed pulse echoes which fall -ntotha range brn «i 
„h JTa rawTeSctk^fe first reported. Since a radar ol the type considered might emptay up to 2000 range tons per 

!Sa Also, multiple detections not we,, isotated in azimuth wfll. for Doppter process^ purposes, be 

dweXT The initial rtrike. received at time t (0>. will be received al some arbitrary phase angle 9 (0) ™ 
grange -transn^erWenoy^ 

(p) and (q) pans of the compressed video pulse in the range gate m whu* the target echo is reported. 
30 6 = arctan (p/q) (radians). 

However, this measurement is ambiguous: 0 can be reported only in .he range 0 to 2 *. The absolute value of e is 
Unkr S2lder now the effect of staggered PRF. In the example the t1 s.rike echoes span 10 Pflls; hence *e jrvwage 

sequence of ±10% amplitude and with sfrtusoktel distribution is now considered. Typ.ca» real stagger patterns are 

relative phase = F(d) x t(p) x 2n radians 
6CP) = 2it % Mantissa [{F(d) x t(p) + 9(0)}/2ii]. 



35 



AO 



45 



The Dopplor measurement system therefore has available to it: 

expressed video pulses in (p) and (q) channels, stored temporarily and covering one dwell period prior to a 
ss detection being reported. 

A raw detection report giving azimuth or time data Irom which the radar strike sequence covering the target can 
be precisely identified amongst the stored video sequence. 
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A rang© bin report lor Ihe raw detection, so that compressed pulse data falling into this range bin can be selected 
from the stored video. 

Knowledge of Ihe PRY patlem(s) used during the appropriate dwell period. 

s 

Using this data, the Doppler measurement system can identify that a particular sequence of compressed pulse 
echoes represents strikes contributing to the detectba By combining the (p) and (q) parts for each compressed pulse 
into polar coordinates representing ampRude and phase, measured (but ambiguous) values of 6 (p) are obtained 
between 9 and 2*. 

to Referring to Figure 3 the horizontal (time) axis covers the appropriate dwell period and is additionally marked to 
show a time sequence of strikes at staggered PRI, taking the reference phase as 9 (0) = 0, for convenience. (For a 
rotating-beam radar the horizontal axis would also represent beam azimuth and could, if desired, be marked to show 
the beam heeding at which each strike was obtained.) 

The vertical axis represents echo phase. The measured phase 9(p) for each strike lies in the range 0 < 6 < 2*. 

»s These measurements are marked by crosses on Figure 3. 

Each of these phase measurements is of course ambiguous in multiples of ±2Nn, for reasons already explained. 
The unknown Integer N will in general take a different value tor each strike. The echo phases corresponding to alt 
plausible values of N Is marked on Figure 3, for each strike, and these are shown as dots, initially, the position of line 
1 , Is not known but it must start from the origin, which In Figure 3 is taken as [8(0) =0; t<6) = 0). 

20 Provided the available phase measurements are accurate, line 1 should pass exactly through one dot per strike, 

representing the true, unambiguous, relative phase shift. The correct position of this line is shown in Ffgu re 3. Its slope 
(4£0w in 10 mi Hisses) represents +1000 Hi which, as already known corresponds to the example 150 m/sec, i.e. the 
true target approach velocity. 

If it Is attempted to draw line 1 through the origin at any other slope, it will ■miss* dots for at least some strikes 

2S (lines 1 a, 1 b) although there remains the possiba&y of a fates result if ft passes suff icientty close to them, (line 1 c) t The 
probability of this type of error increases if the positions of dots are corrupted by noise or interference affecting the 
phase measurements for each strike. However, the probability of error is decreased by ensuring that the stagger pattern 
is not too regular and by not processing very weak strikes (i.e. those with very low amplitude A(p) measurements) 
whose "dot* positions may be seriously corrupt. Subject to this proviso, &s many strikes as possible should be processed 

30 in order to minimise residual ambiguities. 

The slope represents the approach (ie. racSal) velocity of the target and it will only remain constant if target velocity 
remains steady during the dwell time. The maximum plausible acceleration for a target of interest would typically be 
10G, where G (= 9.81 ms* 2 ) is the acceleration due to gravity at the Earth's surface. During a period of 30 ms the 
velocity of a target subject to 10G acceleration would change by: 

35 

ov = 10 x 9.81 X 0.03 = 0.33 m/sec. 

This is very much less than the nominal Doppler resolution (5 m/sec: see above) of the example radar when using 30 
ms dwell HmeFThus whereas in theory target acceleration would cause line 1 to be slightly curved, the effects of this 
over a 30 ms period will he negligible for all targets likely to be of interest 

Line 1 can thus be assumed to be straight and therefore one empirical method of finding the true position of line 
1 would be to pivot a straight thread or edge at the origin of Figure 3 and to rotate It slowly while observing whether it 
passes exactly through one dot for each strike. An equivalent procedure is relatively easy to Implement by using iterative 
as digital processing to calculate an error function based on the sum of -miss" distances between the thread" and the 
nearest dots, while the thread - is progressively rotated through an angular range representing all plausible Dapple* 
shifts. The "correct" slope win be Indicated by a null in the error function. Successful determination of the slope of line 
1 by this means depends on the following assumptions: 

so a) The pivotal origin is taken at or close to the centre of the target strike sequence for this reduces the effect of 

any residual curvature due to target acceleration; 

b) "Good" echo signals for the whole strike sequence can be obtained from a single range bin. By selecting the 
signals to be processed from the range bin containing detection threshold crossings obtained rear the centre of 

6S the strike sequence (Le. at the central origin) the residual loss dun to target motion in fange during the sequence 

of strikes to be processed can be reduced: 

c) That the measured phase o! the slrike selected as the pivot*! origin is not seriously corrupted by noise or 
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interference. This implies that the strike echo at the selected origin should have the best possible signaMo-noise 
ratio ( SNR): this in turn implies a high amplitude A(p) value for this tfrike. Where the radar has a rotating antenna 
beam, this is best achieved by selecting as origin a strong echo at the centre of the strike sequence, where the 
beam points (almost) directly at the target azimuth, H the radar has a phased array antenna which does not scan 
s the target during the dwell time, the strongest near-centre strike echo should be chosen as the origin. This illustrates 
one use of echo amplitude data. 

d) That all other strike echoes processed and contributing to the slype error function have sufficiently high SNRs 
to create a distinctive null at the correct point of the error function, for reasons already given. Since Doppler meas- 
10 urement is cued by a raw MTI detection, it can be assumed that at least some of the strike echoes are strong 
enough to cross the detector threshold and are therefore useful for Doppler measurement. In general, a better 
measurement will result if more weight Is allotted to strong echoes than to weak ones. This illustraies another use 
of echo amplitude data. 

is e) That the PRF stagger pattern in use is clearly unambiguous tor all plausible Doppler shifts, having regard to the 
number of good (strong) strike echoes available for processing. A stagger pattern may need to be oaretully chosen 
in order to validate this assumption. The simple, regular pattern illustrated in Figure 2 is not Weal in this respect, 
in practice an existing MTt radar clutter rejection stagger pattern will probably have to be used unaltered when 
Doppler measurement is added Fortunately, an existing MTI stagger pattern carefully chosen in order to eliminate 

so Doppler blindness should fpso feet© be efftaten I at el&ninating Dopi tier ambiguities. 

Theconditione necessary forthe Invention to provide efficient, unambiguous Dopplerexlraction can be. and usually 
will be, met in any efficient type of MTI radar. -« «™itarf to 

Having explained the principle of the invention, one practical implementation will now be desenbed as appijeoto 
ss a typical MTI radar. The Invention Is not limited to the particular implementation described below Those skilled Inthe 
art will be readily able to apply the Invention to a variety of different MTI radars and in alternative ways with respect to 
points of practical detail. , 

With reference to Figure t. the Doppler management unit DMU, « triggered each time the MTI detector unit 6 
reports a raw detection. The raw detection report identifies the range b.n in which the detection was obtained. 
so The detection process a Iso provides information from which the existing plot extractor unit 7 calculates the azimuth 
position of the target. In the case of a radar with a rotating antenna whose beam (of finite width) scans eorossthe 
target during the dwell lime, the position of the target defines the centro of the target scan, ie. W**^™** 
(or azimuth) at which the beam was pointing most directly at the target and may be presumed to have returned I the 
strongest strike echo, is also provided. Thus the reported target position is expected to correspond with the middle of 
ss a strike sequence although'* will not. in general, coincide exactly with the nominal beam * ■^fffJV 

The digitised (p) and (q) video signals are ted continuously intothe a cycle video store (CVS).B. A given MTI radar 
wai be designed to use some finrte maximum dwell time; the size ot IheCVS is sufficient to retain v.deo s.gnals received 
during this dwell period before they are cyclically overwritten, typically about 2 MByte of RAM. 

Receiver blanking during radar transmitter pulses divides the dfcitised video signal into ^'°™«* 3, 22 
49 shortertnan the radar PRI. Range bin numbering is thus inherent* the sequenceol ^^^f^fj^ 
section. The DMU receives PRF-related synchronisation signals from unit 3 of the exists radar and uses hew to 
time-tag the start of each section ot video stored In the CVS. This enables the exact reception time of every sample 
stored in the CVS to be calculated when required. _ hk! „„., 

When the MTI detector 6 reports a raw detection, intonation defining detection range ^ 
4s position is received by the Doppler management unit DMU B. which Identifies the subset of ^ video sections 
presenting one dwell period centred on the reported target position. From these sections the DMU £^ <rt* 
(q) sample pairs whose (time) positions within the sample sequence indicates that they correspond with the range bin 

reported lor the detection. „. . . n „ , h _, thev 

The DMU 6 copies these sample pairs to a temporary processing store accessed by the DEU 10 so that they 
cann^ be cyclically overwritten before they have been processed. Each of the (typtealy 20) sample pairs copied is 
tagged with the start time-tag of the video section from which it has been taken. flmo |„ ud e S 
The oeu 10 converts each (p) end (q) pair of video samples from Cartesan to pofer format, ie. from amplitudes 
cameled in two Quadrature channels into amplitude and phase [A(s), 0(e). for sample (5)]. , 
^SSSZSZ [ selects the origin. This will normally be coincident with the time and phase of the central samp e 
pair from the range copied, ie. corresponding closely with the reported target posdion. H an even number of sampte 
pairs nave been copied the origin coincides with the centra, pair whose strike echo has M «^ 
measured phase (O S 9 < 2a) of the chosen origin sample is taken as the reference phase (6(0,) tor Doppler extraction 

purposes. 
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The OEU 10 selects tho subset of strike samples which It wfli p>oces$. For strong, isolated detections this wQ! 
normally be the lull set of strike samples copied across from the CVS 8 and covering the tuH radar dwell time. However, 
where the rada r has a rotating b earn antenna which scans past the target heading, early and late strikes in the sequence, 
obtained v/a ihe edges of the beam profile, may be too weak (tow measured amplitude A = poor SNR) to contribute 
5 usefully to the quality of Doppler measurement Also, if two separate detections are reported close together and In 
adjacent or overlapping range bins, the Doppter shift of each target taken singly wifl be better measured by processing 
only strike samples on the "outer" sides of the estimated target centres. 

The DEU 10 therefore employs a sample selection algorithm and rejects weak or doubtful strikes, identified ac- 
cording to simple threshold rules, in order to minimise possible corruption of the Doppler measurement process. In 
io general, the maximum possible number of 'good" strike samples should be retained for processing because this in- 
creases the probability of obtaining clear, unambiguous results. 

The time {-) dt ($) (relative to the chosen origin) ol each strike sample (S) retained for processing is calculated by 
subtracting their lime tags from that of the origin strike. Now. for an arbitrary Doppler shift F(d) the expected phase 6 
(s) of each strike sample in the copied sequence can be calculated: 

e (s> = e (r) + 2ji x F(d) x dt (S) radians. 

The expected phase can be reduced to the range 0 s 8 < 2n as follows: 

20 

6 (s') = 2s x Mantissa [{8 (r)}/2* + F(d) x dt (S)] radians 

Ideally, the (reduced) expected phase 6(a*) for every atrfce processed will exactly match the measured phase 9(S) 
25 when, and only when. 0 (e*) has been calculated using a value of F(d) which exactly matches the true, unambiguous 
Doppler shift of the target which has been detected. In practice, due to noise, Interference, and fluantisation (sampling) 
losses, the match will never be exact; however, the error function should show a deep null at a point corresponding to 
the 'correct* value of F (d). 

The error function is essentially the weighted sum of the "miss" distances d(S) between 0(S) and 6(s') as F(d) is 
30 varied, in general: 

d(S) = l0(s < )-e(S)l Hence 0 £d<S)<2* 

35 However, in cases where 6(S) is close either to zero or 2jt a very small degree of measurement corruption could 

add almost 2* the value of d(S) calculated according to the formula given above. In order to deal with such cases a 
modified parameter cT(S) is used instead of d ($}: 
if 

40 T d(S)£ *.d'(&) = d(S) 

if 

d {$) > Jt. d 1 (S) = 2* - d (5) Hence Oid'(S)<it. 

Thus a typical formula for the error function is: 

50 Error Function = 1/n x Sum (W x {d* (S)} *w], 

for each plausible value of F (d) where n is a normalising factor, W is a weight related to sample amplitude and w is a 
weight determining the power law for the miss distance. 
£5 Suitable values for weights W and w must be selected to suit the application. Normally w = 2 is a suitable value 

for mosl applications. The purpose ol W is to reduce the contribution* of weak samples in order to reduce the effects 
of noise and interference carried within such samples. Where the rariar uses a rotating antenna whose beam scans 
in a known manner across the nominal target position the weight W can be a function ol the off-target scan angle, so 
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thai samples in the centre <* toe copied sequence cany more weight than those near the ends ot the sequence. 
However, il Is advantageous to make the weights easy lo calculate so th-ai the complete error (unciion can be generated 

very quickly. _ _ „ 

The OEU evaluates the error function iteratively over the lull range of positive and negative Doppler shirts F (d) 
s plausible for targets of interest For military aircraft radial velocities < V) in the range ±1 000 m/sec may need to be 
covered Where this is the case, since: 



F(d)=2xV/(X) 



JO 



the range of F (d) to be covered would be ±6.7 kHz in the case of the example radar where X = 0.3 m. 

The error function must cover the Doppler range in iterative steps which are smaller (and preferably quite, a lot 
smaller) than the Doppler resolution which, as already explained, is inversely proportional to the dwell tone covered 
by the strike sequence actually processed. The nominal Doppler resolution corresponding to a processed dwell time 
ts of 30 millisecs is 33.3 Hz {as already calculated). The absolute minimum number of steps required to cover a Doppler 
shift range of ±6.7 kHz would therefore be: 

Minimum number of iterations = 2 x 6700/:i3.3 = 400. approx 

20 Rgute 4 showa an example of an error! unction generated in this way (but using more than 400 iterations). Numerous 
secondary nulls, representing residual PRF ambiguities, can be seen in the error function but the unique, deep null 
indicates the true measured Doppler shift F (d) ; from F (d) the radial velocity of the target is easily calculated by the 
formula already given. . , . 

ss In order to control false alarm rate a null depth threshold is set. Any null not reaching this threshold is rejected. 

This means that, for some raw detections, no Doppler measurement will be reported. In the large majority of cases 
thiswillapplywIytovecyweakdetectiansflowSNR).^ 

of this effect is explained in a later Section. _ . .. 

In practice, use of the theoretical minimum number of iterations, as defined above, may not be sufficient to ensure 
ao that a unique, threshold-crossing null's found in the error function. ev*n for some 'genuine' targets. This is especially 
the case where: 

the prf stagger sequence is inefficient (too short, or loo regular) 
relatively few 'good* strike samples are available; or 
ss SNR is relatively poor 

These types of deficiency may cause false or ambiguous nulls to appear in the error function, as in Figure 4. In 

general, the true null should be deeper than any of the false or ambiguous nulls (as in Figure 4) buts.nce the uxienul 

maybe relatively narrow, its truedepth may not be measured if acoarae iteration step is used. Appreciable Inwement 
40 is obtained by reducingthe step so that approximately threetimesthe minimun\m>mberof iterattonsare us^Jn crneai 

cases, further null improvement is seen when usingaround f 0 times the minimum number. Figure 4 uses2.O0O steps 

each equivalent to 1 .0 m/sec in order to obtain a clear-cut result In practice, selection of an optimum step size must 

take account of the processing time required for very large numbers of iterations. 

one useful stratagem is to generate the complete enor function relatively quickly by using a fairly coarse step. 
as threshold it to find coarse nulls (of which there will typically be several) and then to noajcuUe small « ™ 

function In the vicinity of these nulls, using a relatively fine step in o-der to determine which ot the coarse nulls is m 

fact the deepest and whether it crosses the final acceptance threshold. 

A variant of this stratagem once the approximate positions of candidate nulls have been determined .s to, parfomi 

a weighted 'leastsum of squares-fit (using known algorithms) to determine the absolute minimum null value obWinabte 
so tor eih candidate null region. One advantage of this variant * thai the 'beet fit' algorithm lS not b^^. ^y a 5 - 

sumption that the reference phase 6 (r) of the central strike has been measured wuh negligible error. Th*aliows deeper 

and more reliable null values to be found in cases where SNR is relatively poor 

Those skilled in the art will realise that there are Other known methods, not described here, for refining tn>„ type 

of measurement. Such methods may. of course, be applied H desired. In concept, Do»'^ ext^ perton^ce rnay 
ss be improved to the point where it equals or aven exceeds thai of a good Doppler radar having a s.milar antenna and 

fren3er1a^?»n practice, cost and processing speed resets limit the refinements which can usefully be 

^ The DEU 10 reports the Doppler measurement as a supplement to the raw detection report. Since the OMU 9 
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does not trigger Doppier extraction until a raw detection is reported, the Doppier report will lag the initiating detection 
report by the time required to perform the Doppier measurement. Fc example. 1000 error lunction iterations for a 
20-strike sequence requires 20.000 calculations each requiring, say. 50 basic operations performed at 50 MHZ clock 
rate. The complete Doppier measurement sequence would thus typically take 30 rnillisecs (ie. comparable with the 

5 dweD time) so the Doppier measurement system can just cope with a series of detection reports spaced intervals of 
beamwidth (ie. comparable with the nominal azimuth resolution of the radar antenna). 

In order to unite the Doppier measurement with the raw detection r sport to which it refers, relay of the latter to the 
plot extraction unit 7 is delayed (by approximately one dwell period) until Doppier measurement is complete. One way 
to implement this is for the DMU 9 to store the raw detection data temporarily during the Doppier measurement period 

io . so that detection and Doppier data are always passed on together. Those skilled in the art will know that there are 
alternative ways to deal with Doppier measurement latency. The determination of an approach velocity by the DEU is 
used by the plot extractor 1 to verify that the target is of interest. 

In cases where MTI subtraction has enabled the detector to report a provisional target detection within a region 
of strong clutter, application of the invention to the pre-MTI signal may find and report a Doppier measurement relating 

r5 to the (stationary or slow-moving) clutter rather than to a moving target of interest. Such cases may be instantly rec- 
ognised by the zero or very low value of trie velocity reported. In such ca&es it is necessary to perform further processing 
according to the Invention in order to find and measure the Doppier velocity of the wanted moving target 'hidden' in 
the strong clutter echoes. (Application of the invention to post-MTI signals may be ineffective because the subtraction 
process may cause pnase corruption when PRI stagger is used.) 

2° In such cases the invention can be applied as previously described to pre-MTI stored signals. If this fails to find a 

plausible target velocity it wai instead report the mean velocity of the clutter, which clutter may overlay an appropriate 
expected target position, for example when chaff has been deployed. An estimate of the clutter amplitude is made by 
the Doppier processor (1 0). For strong clutter this can typically be taken as 90% of the signal amplitude in the region 
of the expected target. The Doppier extraction processor digitally generates a synthetic clutter echo, reproducing the 

85 pulse-by-pulse phase differences corresponding to the estimated amplitude, the mean clutter velocity, and to the known 
PRI stagger pattern. The correct absolute phase of the clutter signal ie not known at this stage of the processing, so 
the synthetic signal is generated at arbitrary absolute phase. 

The Doppier extraction processor then performs phase-coherent subtraction of the synthetic signal from the real 
stored signal, measuring the amplitude of the residual difference resulting from the subtraction. This subtraction is 

30 repeated iterative ly while the absolute phase of the synthetic signal is altered in steps. The phase step giving the 
smallest residue is selected. Coarse and fine steps may be used to reduce the total processing load with further iterative 
subtraction possible to minimise the residual signal amplitude and thus maximises the cancellation of the clutter signal. 
The Doppier extraction processor then assigns a target approach velocity procesed from the residual containing target 
returns (and system noise) but no longer heavily overlaid by clutter. Provided the SNR of the residual signal is satb- 

35 factory, this will result m correct measurement of the wanted target velocity, free of clutter. 

Those skilled in the art will recognise that many of the above processes could, if desired* be performed using 
known last Fourier transform (FFT) techniques instead of the simple, but slower, techniques outlines in the exemplary 
explanation given above. However, the techniques exemplified are fully effective and have been tested using target 
and clutter echoes received by a real MTI radar under normal operating conditions. 

ao The invention has been deseffoed above with reference to an existing MTI radar However the invention may apply 
equally to other types of radar and protection is sought for alt such radars within the scope of the following claims. 



Claims 

4S 

1. Radar apparatus comprising: 

a transmitter (4) arranged to transmit a series of pulses; 

a receiver (5) arranged to receh/e and store successive echo signals from potential targets; and 
so a detector (6) arranged to identify echoes of interest by means of a thresholding process applied to data 

originating from the received signals characterised in that a Doppier extraction processor <10) is arranged to 
process stored data corresponding to an echo of interest identified by the detector (6) in order to determine 
the approach velocity of a potential target corresponding to that echo. 

ss 2. Apparatus as claimed in claim 1 wherein the detector comprise* means for subtracting successive reoe ived echoes 
and thresholding the difference signal obtained to determine an echo of interest. 

3. Apparatus as claimed in claim 1 or 2 wherein the stored data is identified chronologically. 
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4. Apparatus as claimed in any preceding claim wherein the Doppler extraction processor pO) is arranged to receive 
from the stored data the relative phases of successive received ec hoe*, compare a range ot velocity values with 
the phase data and identify the velocity value that best fits the ptaise data. 

s 5. Apparatus as claimed in claim 4 wherein the range o1 velocity values is limited toa plausible range for a particular 
target type. 

6. Apparatus as claimed in claim 4 or 5 wherein the Dopplor extraction processor (10) is arranged to determine an 
error value for each ol the range of velocity values, the set of such values constituting an error function, and identify 

10 nulls in the error function. 

7. Apparatus as claimed in claim 6 wherein the Doppler extraction processor (10) determines a finer range of velocity 
values tor a section of an error function about a previously identified null. 

is 8. Apparatus as claimed in claim 6 or 7 wherein the Doppler extracts processor (10) assigns an approach velocity 
to a echo of interest when a null in the error function corresponding to that velocity is below a predetermined 
threshold. 

9. Apparatus as claimed in claim 6 or 7 or B wherein the Doppler extraction processor (10) applies a weighted "least 
20 sum of squares Tit" to a null to determine the minimum of the null 

10. Apparatus as claimed in any proceeding claim wherein the transmitter (4) varies the pulse repetition interval of the 
transmitted signal. 

25 11. Apparatus as claimed in any proceeding claim wherein the Doppler extraction processor (10) only processes ech- 
oes which correspond to a range bin of the radar in which the potential target of interest is first identified. 

12. Apparatus as claimed in any preceding claim wherein the output of Doppler extraction processor (10) is used to 
verify the target as being of interest if a non zero velocity can be r letermined. 
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U. Apparatus as claimed in any preceding claim further comprising means for identifying slow moving clutter at least 
in part in dependance on the approach velocity determined by thn Doppler extraction processor, means for esti- 
mating the amplitude of the signal returned by said clutter, means for generating a synthetic clutter signal and 
subtracting this Irom the received data corresponding to said clutter to obtain a residual echo signal, wherein the 
3s Doppler extraction processor determines the approach velocity of a target corresponding to the residual echo 

signal. 

14. Apparatus as claimed in claim 1 3 wherein the synthetic clutter signal is subtracted from the received signal at a 
plurality of different phase shifts In order to determine the relative phase which provides a minimum residual signal, 

40 whereby the Doppler extraction processor determines the approach velocity of a potential target corresponding to 
the minimum residual signal. 

15. A method of determining approach velocity information from a moving ta/get indicator (MTI) radar, the method 
comprising storing successive returned radar signals, retrieving stored data relating to a target oi interest identified 

4* by the MTI radar, and processing that data to determine the approach velocity of the corresponding target. 

16. A method of determining the approach velocity ol a target, the method comprising: 

transmitting a sequence of pulses of electro magnetic radiation; 
bo receiving and storing data from successive echoes from potential targets; 

identifying echoes of interest by applying a thresholding procass to data originating from successive received 
echoes from a potential target; and 

for an echo ot interest determining the approach velocity of that target 
55 17. A method as claimed in claim 16 or 1 6 further comprising: 

retrieving from the stored data the relative phase of successive echoes from e target ol interest; 
comparing a range of velocity values with the phase data; and 
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identifying the velocity value that best fiia the phase data. 

18. A method as claimed in claim 17 further comprising determining an error value tor each of the range of velocity 
values, the set of such values constituting an error function and identifying nulls in ihe error function, 

5 

19. A method as claimed in claim 18 comprising the further step of determining the error function for a finer range of 
velocities values about a previously identified null in the error function. 

20. A method as claimed in claim 18 or 19 comprising assigning an approach velocity to a target of interest when a 
10 null in the error function is below a predetermined threshold. 

21. A method as claimed in any one of claims 18, 19 or 20 further comprising applying a weighted "least squares fit* 
to a null to determine the minimum of the null. 

t5 
20 
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